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3  QUANTUM  TRANSPORT  AND  MONTE  CARLO  TECHNIQUES 


Since  the  end  of  the  previous  contract 
< PAJA4U- 83 -C -0039) ,  a  major  inprrvement  toward  a  Monte  Carlo 
solution  pf  the  Liouville  quantum  transport  equation  lor  the 
density  matrix  y  (.?.£')  of  an  ensemble  of  electrons  in  a 
crystal  interacting  with  phonons  and  with  an  external 
electric  field  has  been  accomplished.  It  consisted  in 
moving  from  \t  finite -difference  formulation  of  the 
tiiae-cvoluLion  equation  of  the  density  matrix  to  an 
.iterative  expansion  of  its  integral  form: 


o  ’  .  0 


where  Hep  is  the  electron-phonon  interaction  hamiltonian,  k 
and  <n<~,  >  are  the  electron  and  phonon  coordinates, 
respectively . 


From  the  analytical  »>Jevelcpment  of  the  various 
perturbative  corrections  H  XL),  it  is  found  that: 


(i)  the  zero-order  term  is  the  ballistic  translation  during 
the  interval  <0,t.)  due  to  the  action  of  the  external 
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electric  field; 

(i.i>  if  we  want  to  evaluate  mean  values  of  electron 
operators  A<k)  that  are  diagonal  on  the  set  of  states  li?  > 
used  to  define  the  density  matrix,  only  the  diagonal  terms 
of  are  involved  in  the  trace  operation: 

=  Ttete-  (2) 


It  can  be  proved tanalvtically  that  every  odd  perturbative 
correction  £  mY£.E'ifcV  is  zero  if  the  initial 

condition  for  the  total  densi ty  matrix  of  the  electron  and 
phonon  system  is  decoupled  into  a  product  of  the'  electron 
and  the  phonon  part. 

(iii)  all  the  even  perturbative  corrections  ^^(t)  can  be 
written  as  sums  of  terms,  and  each  term  can  be  described  as 
a  product  of  real  and  virtual,  phonon  absorptions  and 
era  i ssions . 

The  numerical  procedure  for  the  solution  of  Bq.fi )  is 
being  developed. 

An  oral-  communication  about  this  recent  work  (abstract, 
here  enclosed)  has  been  submitted  and  accepted  at  the  next 
..International  Conference  of  Semiconductors .  It  will  be 
presented  in  Stockholm  this  coming  august. 

A  preliminary  paper  on  this  subject  has  been  presented  and 
discussed  at  the  recent  Workshop  on  the  Physics  of 

U1  trasmai  i  and  Quantum  Structured  Devices  CJ'empe,  decenber 
1985).  The  preprint,  is  also  enclosed  (C.  dacoboni  et  al ,  to 
/  be  published  in  Superlatt.ice  and  Bicrostructures) . 

•  ‘ 

HOT- PHONONS  IN  SKK I  CONDUCTORS  ,  A>  •*->-’  -  " 

c 

He  have  recently  started  a  project  to  study 
non -equilibrium  phonon  effects  in  semiconductors  using  Monte 
Carlo  techniques.  It.  is  well  known  that  under  the 
application  of  a  sufficiently  strong  electric  field,  the 
mean  energy  of  the  carriers  exceed  the  thermal  one.  Uithin 
this  "hot  electron"  regime,  the  phonon  distribution  can  be 
also  shifted  from  its  thermal  equilibrium.  This,  in  turn, 
can  significantly  influence  the  low-temperature  transport  i.: 
semiconductors . 

The  research  is  done  in  collaboration  with  Prof.  Kocevar 
of  the  Physics  Department  of  the  University  of  Graz, 
Austria.  The  initial  stage  of  this  collaboration  has  focused 
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1  on  the  development  of  a  novel  Monte  Carlo  procedure  that 

includes  phonon  disturbances  in  the  analysis  of  steady-state 
high  field  transport  in  p-Ge  at  low  temperatures.  The 
algorithm  is  based  on  the  iterative  solution  of  the  coupled 
Boitziaann  equations  for  carriers  and  phonons.  At  each 
iterative  step  the  phonon  distribution  is  taken  from  an 
analitical  model  parameterized  by  the' carrier  mean  energy 
and  drift  velocity,  and  is  feeded  int.o  the  Monte  Carlo 
simulation  at  the  following  iteration. 

The  first  results  have  been  presented  at  the  4th 
International  Conference  on  Hot  Electrons  in  Semiconductors,  •  , 

Innsbruck  1985  <see  P.  Bordone  et  al .  ,  Physica  134b,  169  i 
<1985)).  At  the  lowest  electric  fields,  the  phonon 
I.  perturbation  produces  an  increase  of  the  drift  velocity 

*  (drag  effect),  uhile  at  the  higher  fields  the  frictional 

effect  of  the.  enhanced  phonon  population  dominates,  leading 
to  a  decrease  in  the  drift  velocity.  A  more  extensive 
discussion  will  be  submitted  for  publication  shortly. 

A  full  Monte  Carlo  technique  for  both  carriers  and  phonons 
is  now  under  development  for  the  case  of  the  relaxation  of 
photo-excited  electrons  in  bulk  GaAs  and  in  GaAs-All’aAs 
quantum  wells.  A  copy  of  the  abstract  submitted  and  accepted 
at  the  next.  International  Conference  on  the  Physics  of 
Semiconductors  (Stockholm,  Aug.  1986)  is  also  included. 
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